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Abstract. Rashba spin-orbit coupling together with electron correlations in the 
metallic interface between SrTi03 and LaA103 can lead to an unusual combination 
of magnetic and orbital ordering. We consider such phenomena in the context of the 
recent observation of anisotropic magnetism. Firstly, we show that Rashba spin-orbit 
coupling can account for the observed magnetic anisotropy, assuming a correlation 
driven (Stoner type) instability toward ferromagnetism. Secondly, we investigate 
nematicity in the form of an orbital imbalance between d xz / d yz orbitals. We find 
an enhanced susceptibility toward nematicity due to the van Hove singularity in the 
low-electron-density regime. In addition, the coupling between in-plane magnetisation 
anisotropy and ncmatic order provides an effective symmetry breaking field in the 
magnetic phase. We estimate this coupling to be substantial in the low-electron-dcnsity 
regime. The resulting orbital ordering can affect magneto transport. 
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1. Introduction 

Rashba spin-orbit-coupling (SOC) effects have been mostly studied in weakly interacting 
systems such as semiconductor heterostructures designed for spintronics applications PQ, 
and occur in two-dimensional systems without mirror symmetry [2]. However, the effects 
of Rashba SOC in a two-dimensional system with strongly interacting electrons, found 
for example in interfaces, are emerging as a new frontier. There is thus a pursuit for new 
emergent phases of matter in this regime, both theoretically[3] and experimentally p[|, 
with the electron gas at the interface between the two non-magnetic insulators LaA103 
and SrTiOs (LAO/STO) the widest-studied such example. The observed ferromagnetic 
instability at this interface [5] could well be of Stoner type, which suggests that electronic 
correlations may be enhanced due to low dimensionality and poor screening at low 
densities. Hence, combined with its demonstrated tunabilityP, [7], the LAO/STO 
interface is an ideal testbed for physics of Rashba SOC in correlated electron systems. 

Recent observation of magnetic anisotropy may signal further richness in the phase 
diagram of the interface. Specifically, Bert et al. [8] and Li et al. [9] observed strong 
in-plane preference for magnetisation [see figure [U (a)]. Bert et al. attributed this 
observed anisotropy to the shape anisotropy of the interface: energetic bias towards a 
certain magnetisation direction e.g., along the longest axis of an ellipsoid, driven by an 
anisotropic demagnetisation field. However, for ultra thin films consisting of only a few 
atomic layers this effect is subdominant next to microscopic effect s[T0j [TT]. 

If the shape anisotropy is negligible in the interfaces, the dominant source of 
magnetic anisotropy would be spin-orbit coupling. However, the typical alignment 
of the spin with the largest angular momentum in itinerant d-electron systems [T2] 
would predict an out-of-plane magnetisation upon occupation of d xz and d yz orbitals, 
in disagreement with experiment s |13 1 [9j [8]. In this work, we show that Rashba spin- 
orbit coupling leads to the unusual circumstance of an anisotropic (spin) susceptibility. 
Assuming Stoner ferromagnetism close to a van Hove singularity near band edges due 
to SOC, we argue that this in turn leads to a magnetisation anisotropy. 

We further investigate the possibility of nematic order in the form of orbital ordering 
between the d xz and d yz orbitals, since the large density of states near the band edge 
also leads to an enhanced tendency towards such ordering. This order can couple to an 
in-plane magnetisation anisotropy and we show here that a consequence of the Rashba 
SOC is a strong such coupling in the low-density regime near the band edge. This 
implies that orbital ordering accompanies the magnetic phase and leads to an additional 
magnetisation anisotropy within the plane. 

2. Anisotropic Susceptibility 

We use a three-band model for the Ti t2 g orbitals in the xy plane to describe the 
electronic structure of the interface |14j. For now, we ignore the atomic spin-orbit 
coupling in order to gain more analytic insight. In the presence of an external magnetic 
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field H, the Hamiltonian reads 

H — Ho + ^^> c ~~ HbH ■ S . (1) 
Here, "Ho is the hopping Hamiltonian 

^0 = X £ik^ C ]ks C iks ( 2 ) 
i,k,s 

with bare dispersions = k^./2m x + k y /2m l y — \x\ and cj ks creates an electron in band 
I = (1,2,3) =(d xz , d yz , d xy ) with momentum k, and spin s. We use mass parameters 
from reference [15] : the light masses m = m x = m y = m x = m y = 0.7m e and the heavy 
masses M = m y = m x = 15m e with m e the electron mass. The chemical potentials 
are related by \i = \L\ = fi 2 = ^3 + A and we use in the following A = 50meV. In 
the Zeeman term [the last term in equation (0Q)], S = J2i,u, s ,s' c lks^ss'C lks , is the total 
spin with a being the Pauli matrices. Finally, 7-if oc is the Rashba SOC at the interface 
due to the absence of the in-plane mirror symmetry, which can phenomenologically be 
introduced as a relativistic effect due to an electric field E in z direction: The spin of an 
electron moving with velocity v couples to an effective magnetic field (E/c x v). Since 
the velocity v; = /9k for an electron in band / is k-dependent, the Rashba SOC is 

^IL = Q E 9lk- (cL^s'W): ( 3 ) 
Lk,s,s' 

where = (vi ;y , — Vi >x , 0) and the overall scale a ~ 10~ n eVmp!6]. The Rashba 
term %^ oc changes the zero-field bandstructure by splitting the spin degeneracy of the 
individual bands, as shown in Fig. W[b). Notice how different bands have different band- 
edge configurations: a ring of lowest energy momenta for the two-dimensional (2D ) d xy 
band and a saddle point for the two quasi-one-dimensional (quasi-lD ) bands stemming 
from the d xz and d yz orbitals. This difference is due to the isotropic (anisotropic) 
momentum dependence of the Rashba coupling g^ for the 2D (quasi-lD ) bands and 
leads to different types of divergences in the density of states p(e) at the respective band 
edges: A 1/ y/e divergence at the bottom of the d xy band and a logarithmic divergence 
near the bottom of the d xz , d yz bands. In the regime of low electron densities, this system 
can thus have instabilities to broken-symmetry phases even for weak interactions. 

Now we turn to the impact of 7^ c on the in-field (H 7^ 0) bandstructure, which 
leads to one of our main results: the anisotropy in the bare uniform susceptibility 
near band edges. In order to calculate the bare susceptibility, we first diagonalize the 
Hamiltonian ([T]) for each momentum k to obtain the in-field spectrum £\t U (H), where 

v — 1, ... ,6. Then, the (diagonal) susceptibility is given through \% = 

u = Q/N the grand potential per lattice site, 

iV^ 1 4Tcosh[^(#)/(2T)]^ dHi J {m)2 \)\ S= » W 

where T is the temperature and np(£) is the Fermi distribution function. In the absence 
of the Rashba term, the Hamiltonian ([1]) is diagonal for the spin-quantisation direction 
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Figure 1. (a) The distribution of the dipole- moment direction in terms of the 
angle from the z axis as observed in reference [8] (reproduced with permission), (b) 
Dispersion e(k) and density of states p(e) of the three-band model equation ([]} along 
k x . The two quasi-lD bands are shifted by A = 50meV compared to the 2D band. 
The inset shows a spin-texture along Rashba-spilt parabolic bands. 



parallel to H and ^(H) is linear in H. Hence, only the first term in equation (j4j) 
contributes to the susceptibility: the usual Pauli susceptibility. Since the Rashba term 
H^ oc does not commute with the Zeeman term, £ kl ,(i?) is not linear in H once T-if oc is 
present and the second term of equation (0J, the so-called van Vleck susceptibility, 
becomes non-zero. The direction-dependent balance between the two contributions 
determines a possible anisotropy Xz ^ Xx- 

As the Hamiltonian in equation ([1]) is block diagonal in the orbital basis, the total 
bare susceptibility is the sum of contributions from each orbital I. We start with 
the contribution from the quasi- ID orbitals. The two d xz bands have dispersions: 



2m 



7.2 
K v 

1M 



H ± \{ag xz ^ - fi B H)\. 



They contribute to the total susceptibility through 

1 



and 



xf(xz) 



xf(xz) 



k,± 



4Tcosh[£^(0)/(2T)P 



/4£[i 

k 



\dxz,k\ 



(5) 



(6) 



(7) 



with g xz y the unit vector along g xz ^. For T — > 0, we substitute k y = (M/m)k y and 
change the sums in equations ([6]) and ([7]) into (cylindrical) integrals. For /i > 0, we 
obtain 



and 



t4M 
2tt 2 



(9. 



XZ,(f>) 



COS z 



+ M siri 2 , 

m 



(8) 



1 - (51. 




Figure 2. In- plane (solid line) and out-of-plane (dashed line) total spin susceptibility 
for the three-band model (a) without and (b) with atomic SOC (a at = 5meV) with 
the gray area denoting the region with anisotropic susceptibility. The scale in these 
plots is given by Xq = l J -B m / 7T an< ^ Xq 11 = Xo V + 2/igVMm/7r (right vertical axes). In 
(a), the contributions of the two quasi-lD bands for Xx are also shown separately. The 
arrows in the insets denote the preferred magnetisation direction, where in (b) also the 
anisotropy outside the 'Rashba' regime due to atomic SOC is shown. 



with g l xz ^ = (sin 0, — cos0, 0) and is the angle relative to the crystalline x-axis. 
Clearly, the total d xz contribution to the susceptibility Xi( xz ) = l^vrnM /it is 
independent of the field direction % or the Rashba SOC strength a. 

On the other hand, near the band edge of the one-dimensional bands, i.e., 
—a 2 /2m < fi < 0, the susceptibilities in equations ([6]) and (GO) yield 

xfM = 4¥ / d<P }f X *M- 2 J l+2/im(cos 2 0+^sin 2 0)/ a r i/2 (lO) 
2tt z J cos z o + — sin d> m 

and 

xfM = / d<P \ ~ { f X ^ ? 2 7 [l+2^(cos 2 0+^ sin 2 0)/a 2 ] 1/2 -(H) 

27r^ J cos + — sin m 

The total contribution to the susceptibility is thus anisotropic at the band edge. 

For the total bare susceptibility, we also need to consider the 2D orbital d xy 
contribution. We can read off Xi( x v) from the quasi-lD contribution Xi( xz ) discussed 
above by setting M = m and shifting the band edge by A. Therefore, Xi( x u) = Xo V is 
isotropic for \i > —A. When contributions from all the orbitals are combined, the total 
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Figure 3. Coupling constant A cx Xxxrj for coupling 77 = n\ — r&2 to — as a 
function of chemical potential for different atomic SOC parameters. Note that A for 
atomic spin-orbit coupling only would in this plot not significantly deviate from the 
zero line (thin dotted line). 



susceptibility shows two regions in the chemical potential where Xx = Xy > Xz- near the 
band edge of the 2D bands and that of the quasi- ID bands. Fig. |2(a) shows the total 
susceptibility for the three bands near the quasi-lD band edge as a function of chemical 
potential, where we shaded the anisotropic region. 

Fig. Efb) summarises two ways atomic SOC impacts the phase diagram: (i) it 
widens the region of anisotropic susceptibility (grey region), (ii) it causes anisotropy in 
the spin direction by aligning the spin with the largest angular momentum when the 
susceptibility is isotropic (white region). The first effect results from adding the atomic 



to the Hamiltonian ([I]) and evaluating the susceptibility fl4]) numerically. For the second 
effect, the atomic SOC is treated as a perturbation in the magnetic phase[12]. This 
again predicts in-plane magnetisation below the grey region where only the 2D band is 
occupied, but a switch to out-of-plane magnetisation above this region. 

In order to compare our phase diagram with experiments, we need to translate 
the chemical potential to a gate voltage. While such a translation is non-trivial, 
Hall measurements under gate-voltage sweeps can offer hints as to where the as-grown 
samples lie. Joshua et al. [13] showed that the interface acquires heavy carriers when 
under a gate voltage, which indicates that the as-grown samples are near the quasi- ID 
band edge. Though the grey region is narrow in Fig.[2](b), the density changes by a factor 
of ~ 2 in this range. This pushes the density upper bound for in-plane magnetisation 
substantially, consistent with experiments [17] . 



SOC, 




^ e lmn ^ C | 



£ks c mks' u ss' 



(12) 
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3. Orbital Order and Magnetisation 

In the regime of low electron densities, the van Hove singularity can promote broken 
symmetry phases in the presence of suitable interactions. While we focused so far on 
a magnetic instability, which could for instance be driven by repulsive intra-orbital 
interactions, we now turn to orbital-ordering possibilities. Even though the three-fold 
degeneracy of the t2 9 orbitals is already broken due to the interface symmetry, the d xz 
and d yz orbitals remain degenerate. A spontaneous orbital symmetry breaking described 
by a non-vanishing order parameter 77 = ri\ — n 2 , with 774/2 the occupation of the d xz /d yz 
orbital, can be driven by inter-orbital repulsive interactions [18]. We first analyse the 
tendency toward such an instability by studying the (bare) nematic susceptibility. For 
this purpose, we introduce the field H v conjugate to 77, which enters the Hamiltonian ([1]) 
through fii/2 = fi ± H n , i.e. it acts through an opposite shift of the chemical potential 
for the two orbitals. The nematic susceptibility then yields 

_ d 2 co d 2 cu 



d 2 cu 



and is given by the contribution of the d xz and d yz orbitals to the density of states. 
Hence, the van Hove singularity near the band edge of the quasi-lD bandsfsee FigfjTJb)] 
allows for a nematic order for sufficiently strong inter-orbital interaction [18]. 

Next, we consider the coexistence of nematic and magnetic order. While a system 
with C4 symmetry does not allow for coupling of rj and \M\ directly, there is an allowed 
tri- linear coupling between 77 and an in-plane anisotropy M 2 — M 2 , as both acquire a 
factor of —1 under C4 rotation. Specifically, a coupling of the form \(M 2 — M 2 )r] enters 
the free energy. Given the Hamiltonian ([1]), we can explicitly calculate the coupling 
constant A. It is given by the generalised susceptibility 

X(XXvxx = (dH v dH 2 \i}\=o' (14) 

which measures the change in the spin susceptibility Xx upon shifting the chemical 
potential of the two quasi-lD bands against each other. Figure [3] shows the result 
with and without atomic SOC in the presence of Rashba SOC. (We find the coupling 
to be negligible without Rashba SOC.) The coupling becomes substantial in the above- 
identified density range with the in-plane preference in the bare spin susceptibility. This 
is due to the inequivalence between the two quasi- ID band contributions to Xx shown as 
dotted and dash/dotted lines in Fig. EJ^a), more specifically the difference in the slope of 
Xx(xz) and Xx{yz). The sign of A determines the relative sign between 77 and (M 2 — M 2 ) 
and whether the majority orbitals will be along or perpendicular to the magnetisation 
axis. Notice how our result shown in Fig. [3] predicts a change of the sign of A over the 
chemical potential range exhibiting the magnetisation anisotropy. 

Given the observed magnetism, we finally investigate the effects of coupling between 
magnetic order and nematic fluctuation associated with nearby nematic phase. In 
particular, we may ask on the one hand how the proximity to a nematic instability 
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influences the in-plane magnetisation, and on the other hand, how a magnetisation 
effects the orbital ordering. Assuming an XY ferromagnet in the absence of a coupling 
to nematic order, the Landau free energy becomes 

f(M x , M y , rj) = f + ^p|M| 2 + b -\M\* + X(M 2 X - M 2 )rj + + • • • , (15) 

where M x (M y ) is the magnetisation along the crystalline x (y) axis. Since we assume 
no independent instability towards an orbital-ordered nematic, a v > C0- Integrating out 
rj by minimising the free energy with 

V = -±(M 2 X -M 2 ), (16) 

the free energy for the magnetisation becomes 

f(M x , M y ) = h + ^\M\ 2 + h -\M\'- ^-(M 2 X - Ml) 2 . (17) 

For finite A, the coupling to orbital order therefore locks the magnetisation along one of 
the crystal axes, leading to an additional in-plane anisotropy. 

For the orbital ordering, the magnetisation anisotropy acts as a driving field. For 
a v very small, i.e., the system close to an instability, we should include additional terms 
to the free energy for t], 

f(rj; M) = |V + ^ + ^ + X(M 2 X - M 2 ) V + ■■■. (18) 

For b v < 0, the system undergoes either a metanematic crossover, or a first-order 
transition at a critical magnetisation. This is in analogy to metamagnetic transitions 
observed in systems close to ferromagnetism[T5] and could here indirectly be driven by 
an applied magnetic field. 

4. Concluding Remarks 

We have shown that the combination of Rashba SOC and atomic SOC leads to an 
electron-density dependent magnetisation anisotropy in LAO/STO interfaces. While 
experiments so far appear to lie in the in-plane-magnetisation region, we predict a switch 
to out-of-plane magnetisation at sufficiently high gate voltages. We have identified a 
regime near the band edge with anisotropic susceptibility as a non-trivial effect of the 
Rashba SOC in a low carrier density system. The high density of states near the 
band edge in principle also allows for a spontaneous orbital ordering and we predict in 
this regime an enhanced coupling between the magnetisation direction and this kind of 
nematic order. This coupling locks the in-plane magnetisation direction to be along one 
of the crystal axes and promotes Ising nematicity. 

Next, we comment on the issue of heterogeneity detected in reference [S]. 
The observed heterogeneity is likely driven by both extrinsic and intrinsic effects. 

X Even if strong correlation drives orbital order there will be no qualitative change in the effect of the 
coupling A in driving the in-plane magnetisation anisotropy. 
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For instance, a recent study showed that strong Rashba SOC can promote phase 
separation [20]. The proposed magnetisation-nematicity coupling has important 
consequences in both extrinsic and intrinsic fronts. On the one hand, oxygen vacancies 
and other spatial inhomogeneities act as a random field for the Ising nematic and in 
turn cause a distribution of moment directions. On the other hand, the reduction 
of the magnetic order parameter symmetry due to the coupling changes the type of 
magnetic textures and their energetics. Furthermore, we expect the coupling to cause 
non-trivial in-plane anisotropy in the magneto transport [§|. Moreover, the sign change in 
the coupling A could be observed through the rotation in the dominant direction upon 
gate voltage sweep. Finally, we note that the range in the density with magnetisation- 
nematicity near the band edge of the quasi-one- dimensional bands has also been shown 
to exhibit critical scaling in recent Hall measurement |13] . 
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